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A^.'M'RACr 

A ■no<L->i in 'aPi i c‘i thsrinnl t?f foots are simulated through us«? 
of a nultilieam plasma li st r ihu t if'>n function is developed and 
i nv<'S t i ga f I' 1 to sr>e» if sol ijt ions which take an initially uniform 
magnet iJ^O'i oiasma to a new uniform r>tate with a different field 
r>r i ent at i on are oonsible. 'I’he momentum cr^nserva t ion integrals 
are ‘"ound to admit two claene'^ of such solutions, but only one 
class exhibit a aporooriate asymptot ic behavior. Extensive numer- 
ical integrations have f.iilel to demonstrate the existence of the 
les i re,'' ;',at ijti ons . 


f 


I 
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1.3 TN’PRr)DUCTIOM 

At present there exists no theory which flescribes the 
structure oT larqe nmplituie plasma disturbances in which the 
magnetic Ciei i is rocatel while plasma flow remains parallel to 
the field both ^ar upstream and far downstream from the 
disturbance. Jumo conditions approoriate to such structures are 


known from MMO 

theory , 

and 

include 

true rotational 

discontinuities in 

which only 

the 

field 

direction changes. 

intermeliare v/aves 

in which 

plasma 

f low 

goes from super to 


sub-alfvonic and which aopear nonevolutionary in MHD (Jeffrey and 
■Taniuti, ld64) , and mixed structures which occur in anisotropic 
olasmas and which exhibit characteristics of !>oth shocks and 
•'otational discontinuities (Hudson, Ju70). 

■'lork with larq-^* amplitude disturbances having the desired 
upstream orooerties has hr>pn carried out by Montgomery (ly59) , 
Saffman (Idhl), Kel log (1304), and Kakutani (1960), and special 
enses are summarizeil loy Tidman and Krall (1971). These authors 
used similar cold plasma models, and found solutions representing 

infinite wave trains as well as solitary waves in which plasma 

* 

an 1 magnetic field parameters ar#» identical on cither side of the 
disturbance. Tn no cases were field rotating solitary waves 
found, and it is relatively easv to show that such waves are not 
r»oscible in col 1 plasma. The closest solutions to those desired 
are Raffman's "quas i sliocks " , in which an initially uniform plasma 
unlerqoe.s e disturbance after which plasma and field parameters 
v/ander ergod icr d 1 v ahiout certain mean values, but never approach 
the uaiform conditions sought here. 
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Tntorojjt in Tiel.i rotating wavon has boon rekin'-lled by the 
availability of nagnetic fielrl lata from spacecraft crossing both 
the terrestrial an<l Jovian magn«‘?topauses . During times when the 
so^ar tie^d is approximately antiparallel to the planetary field, 
the magnetopause r^'sembles a layer through which plasma may flow 
and across which the transverse magnetic field is rotated. Many 
magnetic fielrl signatures are remarkably laminar and resemble in 
polarization and overall, shape the solitons of Saffman and 
others, if somehow only half the soliton solution could be 
isolated. ’’he available data and their relation to current 
theory have heon summarised by Sonnerup (1977). 

The present paper pursues a suggestion by Sonnerup that 
thermal erfects which niigiit oermit field rotating layers could be 
i ntroJuced into the tbrn^ry bv assuming a multibeam plasma 
distribution runction. This work shows that the momentum 
conservation conditions, which do not generally admit field 
rotating solutions in the coll plasma case, do admit two classes 
of such soluti<3ns when thermal effects are inclutled. Asymptotic 
analysis shows that liniform solutions appropriate to the 
downstream conditions are capable of growth for one class *ut not 
the other. 


2.0 nOUATTOMd OP r'HK MUbTIIiPAM TMKORY 

riie analysis is performed ^or a one dimensional case in 
which all variables are functions of z only. Assume the plasma 
tf) consist laf M beams of ions and a single beam of electrons, 
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with tho q'jas i nout ral i ty oon<iition satisfied. The electrons are 
assumed to Follow the magnetic field, while the electric force on 
the ions »s neglected. Although these latter two assumptions can 
be relaxed, they In not appear to affect qualitatively the 
results obtained the waves being sought here. Finally, there 
is assumed to be no current in the z direction. The equations 
describing the svstem are then: 
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( 6 ) 




M< . \J. 1 / 
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^ X f<: u (S 


I 


(7) 



Vi.,. V 






cur^r?m' .^ohic ? 


( 8 ) 


V cintl n an? the electron velocity and densitVf ^ and the 
velocity and donn itv of the ion beatn^ V the ion mass, e the 

clopnentary charqn, n the magnetic field, and c t!ie speed of 


li<]ht. Ouan.-itiee -nt!? the subscript „ reoresent values at some 

o 

vchosen initial ooint. All sums run from I to N, whore N is the 
fuimb*?r of beams. Ttiose equations are obvious generalizations of 
t ie usual two riuiil m<>lol. However, they can i>o shown to follow 
rigo>'ous1y from the VI asov-'Mxwel 1 equations in the case of a 
multibea;n distribution function when looping orbits (ie, negative 
values of V^) are not permitted (Jones 1977). 
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TntroHurinq equations (7) nnrl (8) into equation 


(3) gives 





(9) 


On the other han.l, the x component of the ion motion equation may 
he V 1 .V i tten 



run 1- i m y inq th i r. equation bv and summing over k gives 










'■'omnarinq equatir>MS 


(3) and (10) allows an integral of the system 


N\ 


ok 


4 Ok. ^/( 


4^ 

V7T 




to bo written: 
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where If, const-vnt of the motion. An irientical procedure 
qiives an integral associated with the y component; 


M Z ^tn\< ^Y/f 



( 12 ) 


The z component is handled in a similar but not identical manner » 
y i e 1 U ng : 


h 


r>L 


1/ 

^ojf. 


l/.„ + 


'^K 


1 

s-;r 



(13) 


Equations (11) - (13) state that the flux of momentum is 

constant . 

The ion motion equation (6) and the momentum equations (11) 

- (13) constitute 3 n+ 3 equations in the 3N+2 unknowns B , D , and 

X y 

One ef^juation is redun lant and can be used for checking 


the V^'s. 
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°0?'' POO"« QOAUTV 

nurneric.il results. In practice, equation (13) is best suited to 
this purpose. 

because the ions Peel only the magnetic force, their 

velocity m.ignitu<les remain constant. Thus N additional integrals 

2 2 

can he formed expressing constancy of the terms xk yk 

2 

V Und^r the no loop condition there is no ambiguity in the 

sign of V , , so these integrals can be used to eliminate V . 

/u K Z K 

a 1 together . 


3.0 riNlPORM FTl^LP ‘30LUTT0N 

Th<^ largo jinplitiido waves hoing sought woulrl connect two 

regions in which the flow and Field are asymptotically uniform 

ind aligned with each other. In the single beam case (N=l) the 

only wa'r in which such uniformity can occur is with the single 

velocitv vector parallel to the magnetic field. Then the 

derivatives in equation (5) vanish and the solution is truly 

uniform. Tn this single beam case it is easily seen that the 

momentum equations admit no large changes in field orientatign 

except for the snocial case P = P =0, where an exact field 

X y 

reversal aopears possible. In this case, however, all 
derivatives v.inish identically and there can be no wave. 

In the multi beam system, however, a situation in which 
individual beam velocities are not parallel to the field may 
nevertheless produce a net flow along a uniform field. In this 
case the velocity components perpendicular to the field represent 
a sort of quasirandom d i st r i l)ut ion which simulates the effect of 
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a nonzero Ion temperaturet To quantify this situation, let all N 
beams have the same velocity magnitude, and let each velocity 
vector make the same angle with the field. Call the tangent of 
this angle X , so that is the ratio of perpendicular to 
parallel velocity comoonents. is thus a measure of the 

"temperature" . 

Wlien N=I a nonzero "if violates the uniformity condition, and 
in fact gives rise to the infinite wavetrains of Saffman (1961), 
When N>i, it is reasonable to look for solutions in which the N 
beams carrv the same particle flux, have their velocity vectors 
snaced as regularly as possible about the field, and rotate about 
th«? field at the larmor frequency. With call beams carrying the 
same flux, the quantity n , V . appearing in equations (11) 

(13) is independent of k, so that in a uniform field region those 
equations take the forms; 

^ ^XK ' Conjt/inT 

(14) 




( 15 ) 
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^ - con?T/«'^T . 


(16) 


Furthermore, with a uniform fieli, the z component of equation 
(6) can he multiplier! by summer! over k, and integrated to 
give: 



- cohsta^jt. 


(17) 


In the two beam case equations (16) and (17) constitute two 
equations in two unknowns, whLcli require that » '^z2* 

can only happen in a uniform field if the velocities lie exactly 
Parallel to the field, so the uniform field two beam solution 
collapses to the one beam cold plasma case. With N>2 this 
constraint is lifted, so that N=3 is the simplest case in which a 
nontrivial uniform field solution can exist. 

r'onr. ider the particles of each beam to follow helical 
trajectories lifferinq onlv in ohase. It is necessary to write 
equations "or th'^sp trajectories as functions of z when the 
particles spiral about a field which does not lie along the z 
axis. The situation is shown in fiaure 1. Only a single trajec- 
torv is shown, which passes through the origin. However, all 
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snace is ♦'ille'^1 with traiector ies of particles of this beam, and 
the trajectories are all in phase at a given 7, not at a given 
distance from a plane peroendicnlar to the field. This situation 
is necessary liecause of the condition that the system be one 
dimensional, with all variables functions of z only. Any wave- 
fronts existing in nonuniform field regions lie perpendicular to 
the z axis. 

The trajectories can be described easily in a coordinate 
system x", y, in v;hich z" lies along Q. In this system the 
trajectory of figure 1 is; 

X}/. k^siN(oj%- 



y 
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Of. fj’UXiiiJ 

when* « McVj^ -'♦'H in tho 1 armor ra'Uus, CO* cR/Mc tho larmor 
froquoncy, an-1 V^y the velocity comixjnents perpendicular 

and Piirallol to H, in<1 n phnsr* .inqle For the beam, is 

a time p.ir .meter proporti«ma1 to the position of the guiding 
center alonq the fiel<l. Recausi? the larmor orbits are not 
perpendicu! ir to the t .ixin, V.^l^ changes with orbital phase* and 
t^ Ik not a .simple fiinct’on of 7 ,. Performimj .t rotation gives 
the parameter i zivl tr.^lectory in the x» y, z system: 

- ^i<) 

( 21 ) 

N/. ' (] Q SIN 

A <9 ^!n(co a ' a ) ^ 

wh*-»e ^ I-; tii,» .\>>ple lit'tweef^ the field and the 7. axis. 


n i f ( *' r en t I a t i s.t -Mui.^tions ( .! I ) 


(21) gives the velocities: 
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whoro, iq.iin, ^ = Vf /'^fj • Although these velocities were 

ierivotl for oArtlcles wliose guiiiing center passefl through the 
origin, they v'oulri rli^ncr Lhe the entire velocity fielrl as a func- 
tii^n of 7. if the relatiiin between tj^ and z were known. 

’=;qu.itit)n (?3) nrovides the lesirefl relation, but is 

un for tunat>d y transcendental. Analytic progress can be made by 
=»ssumin ) jG < < I an*' *>xpan ling in powers of 'b . Because the 


y 2 • 

, i t 1 


is 


necessary t*'> carry t*‘‘rms of this tirdr^r if results differing from 
th<^ col l olnstm theory are to be ol>tained. Since terras involving 
t^ in etiuations (21) - (2b) are already multiplied by V » it is 
only necessary to solve equ.ation (22) to first order in r. TO 
df> so, define a time parameter t - cos^^which is strictly 

proportional to 7 am) write ty, = t At,^, where A is 

assumed ‘-.hall crrmoire to one. T)ien noting that Rj - , 
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and letting “ tan 6) , equation (?3) becomes 

Co - Yfi SM t Az-a'tcJ 

Rimplifvinq an 1 •'et a i n i ng only t^^rms of order t this is 
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(27) 

* 


Insert inq ev^uation {^'7) into equations (24) 


velocity solutions to sec<^nd order 



(26) gives the 



cos&l I Y/> ^/3 (j6) 
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Thf*se are explicit functions of z because t - ?/ % as 6) . 

With uniform phasinq = 2 TTV/n) , it is straightforward to 

2 2 

verify that, to terms of order Y , these equations constitute 
a uniform field solution of equations (6), (11) - (13). The 
inclusion of ^ in the smallness parameter means the field angle 
cannot approach too closely to 71/2. In the process, the follow- 
ing useful results emerge: 


Ui c»s e ^ 


( 33 ) 
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whore total initial ion density. 

Before prnceeiiinq, it is interesting to note that this solu- 
tion is not exactly uniforn. There are small wiggles in the 
rielfl, hilt in the three honm theory they are already of higher 
order than N is increased, comparison of equations (28) 
througli (30) \.ith a numeric.al solution of the transcendental 
equation (23) shows that nonun i form it ies in the field drop off 
rapi lly with increasing N. Also, the field vectors are not 


equally and riqitUv st^aced around the field, but wiggle back and 
forth in relatif^n to one another-, 'this is necessary to keep the 


fluv con?', tant and along J3 as 
dons i t ios , vary. 


individual 


hence 


4.0 FI!’ r.n i<o TAT I ON r’ON n i •!■ i ons 

A tMtt iou'ar uniform field solution is <letermine'l completely 
')V the parameters (9 . V// , V , n„, and n^. Of these n^, being 
an initial ('onditifin, and .xre strictly constant even through a 

tuMiiin i I >>r n Pho other three might vary throughout , a 

r'.onnn i form regii’it^ and emerge changed but again constant to 
'haracteii/e a 'lifferent uniform field solution. This possibil- 
ity i ,s joverni'd in part by the equations: 


--P 


VTT 
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M f _- p, 


S7T 


(13) 




(34) 
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where V ir. the velocity maqnitu'le, which is strictly constant. 

F.q'jations (11) (13) are the x and z momentum conservation 

equations and equation (34) expresses conservation of energy for 

the particles interactinq only with the magnetic field. The y 

momentum equation is not inclnde(^l because it has already 

constrained any nev/ uniform Fields to lie in the x-z plane. 

Consider a uniform field solution characterized by the 

parameters ^ o' o' Yo- When n^ and are also specified, 

the constants in equations (11), (13), and (34) are determined. 

Are there then othor sets ( &. ) , representing other 

^3 r 

uniform field solutions, which also satisfy equations (11), (13), 
and (34) with the constants determined from the initial uniform 


field solution? To find such sets, the sums which appear in 
equations (11) and (13) can he evaluated using equations (31) - 
(33). Inserting these latter expressions into equations (11) and 
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(13), in'l «*vi 1 uating tho conntiants and from the initial 
parameter net results in the equations: 



h *io ^,/ caid f !■' - y/s^) -f A 




§L 

Sir 


Uning efjuation ( 34 ) to eliminate , there result two equations 

in the two un!<nownG Y >6 (recall that tan ^ ) : 


M// 1 . v' 4 ), A-/ rA^.Ji-y)- ^ 

ilT 

' U Vo; 

^IT 


( 35 ) 


fit 

-- /7/v), vV^r'6) ( 


C^) 
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Those equ;^tions can b(* maio •'I i mens ion less by defining a 
siirtablo ratio of kinetic to magnetic energy densities. Normally 
this would he the AlFven number, given by: 




W/ti Ml/ 
/S' 



(37) 


!>ut here 


it apnoars that a modified A1 Fv^n number. 


given by: 



0 ~ I ^ ) 



(38) 


Dlays a more fundamental role. When this number is close to 
unity the definition (38) can be shown to be a special case of 
Hudson'r. (l‘)70) modified Alfven number appropriate to anisotropic 
plasmas, and under these conditions will be referred to as the 
Hudson numljer, m^j. The plasma un ler consideration is indeed 
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."\nisotropir; , since the pressure, as determined by the quasi random 
perpendicular velocity components of the N beams, is entirely 
perpen<Ucular to the magnetic field. 

Equation (3 j) can be solved for , and by using the 
modified Airven number t'le two equations become: 


yVrVl 




SIN 

SIN^ 



S/N60 


(39) 



Equation (3v) can be inserted into equation (40) to give a single 
equation for 9. Thin can be solved numerically to reveal two 

Classen of solutions, one in which the field undergoes a near 

t 

reversal ami the other in which the field angle changes only 

nlightlv. There is, of course, also the trivial solution 0 = 

which corresponds to the soliton solutions of Saffman and others. 

More iindern tan ling of the solutions can be obtained by 

aof)roaching equations (39) and (40) analytically. Because the 

one beam theory admits an exact reversal (® = - 0^) when m^ » 1 

and m^ arc indistinguishable with one beam, since V “ 0) » it 

2 

is reasonahln to \nn'< for solutions with m^^ 


close to one and 
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which nearly reverse the fielrl. Let ^ “ tn„^ - 1 an«i 0 = - (0 

p o 

Putting these forms into equations (39) and (40)., and 
performing considerable algebra while retaining only fitst order 
terms in 6 and ^ results in the near reversal solution; 



(41) 




(42) 


Thus the extra degree of freedom afforded by nonzero temperature 
has allowed field rotating solutions to satisfy the momentum 
conservation integrals for values of the Hudson number other than 
unity. The solutions are shock-like in that the temperature 

t 

r i fies acron^ tho iniicate^i equation (42) • With the 
initial Hudson number given by = 1 , the Hudson number 
in the reversed ♦'ield region can be evaluated from equations 
(18), (•*!), and (42), and is found to he m^^^ » I - ^ . Thus the 
Mu Ison numlier goes from above unity to below. However, the 
Al^'ven number, given in this case by m. = m + remains 
above unity. It is unclear how to classify these possibilities; 
they might .ler.crihe intermediate waves hut if so their evolution- 
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ary character is unclear# since it is not known whether 
intermediate waves in anisotropic plasma are evolutionary (Hudson 
1970) . 

The second class of solution can he analyzed when the Alfven 

2 

number is close to sec 0. Tn this case numerical work suggests 
only a small fie^d rotation, so a solution of the form Q ■ ^ 

is assumed. The analysis is complicated by the proximity of 
the trivial root (!) but eventually yields the solution, 

good to tP'cms of order 



(43) 



(44) 




2 ' 

where now (^ ~ m /sec"^ 0 -1. Equation (44) 

o o 

situation is possil)le even in the one beam case 
Evaluation of the modified Alfven number 


shows that this 
where y is always 
in the rotated 


field re j ion gives 
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(45) 


It is useful to compare the new Hudson number with the secant of 
the new field angle. Introducing equation (43) into equation 
(45) gives the result 



(46) 


Thus .1 Miidson number initially less than sec^ (ie, <5 < 0) 

2 

becomes greater than see d in the rotated ^"ield region. In the 
analysis leadin*} to both rotated field situations^ it is assumed 


that 


^.O A^YMPTOT TC SOLUTIONS 



t 


In Saffman'r, cold plasma theory it is shown that solitary 
waves can exist only when th»» A1 fven number, as defined in the 
present oaper , satisfies the condition 
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(Saffman 1061) . Small por turbat ions of uniform field solutions 
outside this range are incapable of growth. Both field rotating 
possibilities in the multihoam warm plasma theory take the 
modified Ai fven number from a region of permissible Alfv/n 
numbers in Oaffman's theory to a region where asymptotic growth 
is not possible. However, the Al fv4n number may remain within 
the permissible region. Saffman's results cannot be applied to 
the multi beam theorv i>ecause the nonzero temperature and 
associated pressure anisotropy creates a distinction between the 
Alfvln and modified Al fv(?n numbers. It is necessary to perform 
an asymptotic analysis of the multibeam theory to see whether the 
waves admitted by the momentum integrals can actually exist. 

Tt is convenient to make the equations dimensionless through 
tlie "o1 lov^^i^g definitions: 






L 


Q 




€ fip 
M c \/ti^ 
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Merf* the ^subscript ^ refers to cori'litions in a 
reqion. The ilefinition of b gives imuediately 
is the initial angle between R anf! the z axis 
Icf ini t i on-' , equations (11), (1a), and (6) become: 



oft 








uniform field 
cos 0 where 0 
With these 


( 48 ) 
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Di CtiO 





(4») 


(50) 


wh<'ro X » V. const, int iind thf» ^nnloqous y constant has been set 
to 0 so that the field lies initially in the x-z plane. Let 
^ ok^zn't icquirc this quantity to be i mlependent of k, so 

that Initially riu' bcamr. differ only in phase. Then equations 
(IB) and (4‘)) can l»e si-»lvc<l ft>r thn maqnetic field components to 
a i ve : 




( 51 ) 
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Intro-lucinq thene forms into the x and y components of equation 
(30) gives 



(53) 


aUyj^ 

2~r 






COS 6 (54) 


To sof' whc'th'^r growth from th'' uniform fieli solution is p>ossi- 
hlo, lot y.|^ - V w » whore Vj^ is a uniform field solution. To 
first order in the Wj^'s, equations (33) and (54) become 



d hi u - 

{ -x/t 






/ 


I 


(56) 
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O 

Conservation of onerqv requires that 
constant, so that 


+ u 


Yk 


2 


+ 



remain 





(57) 


Assuninq Y to he small, <1 i mens ionless forms of the uniform 
fiel 1 solution (equations (28) - (30)) can he use<i for the 
Tntrodiicinq these into equations (55) - (57) and using equation 
(57) to eliminate yields, after much algebra, the 2 n 

equat inns 

+ v'/* (n\oi2o(k^ 



» 



( 59 ) 



whore t-ho exponential ronresents either a growth or oscillation 
cormion to all beams, and the 3L|^'s are individual and presumably 
small oscillations at the larmor frequency. When equation (60) 
is introduced into equations (50) and (59) there result 
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The perioilic wiggles represented hy the can be described by 

a four ier series solution of the form 









(63) 


The growth r>r oscillatory nature of the overall solution can be 
determiiu?! l»v introducing equation (63) into equations (61) and 
(j 2), lookirni only at the i nilopondent terms, and summing over 
k, which 
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where 11^ = < !l|^o* ^ne of the fi's in these equations is an 
k 

arbitrary scale factor, so that these are two equations in the 
two unknowns a and the ratio of the B's. Multiplying the first 
equation by 3y^# the second by and adding gives 



Solving equ.it ion (04) for and introducing the result into 

# 

equation (b6) gives 



( 67 ) 
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9 

The root ® 0 is of no interest here because, even if it does 
permit growth, the average values of the field perturbations 
remain zero so that such growth would not permit solutions with 
flifferent field orientations to be connected. The second root, 
given by 


OF 



( 68 ) 


requires that the A'l fven number be close to but above unity. In 
this case the modified Alfven number is the Hudson number and is 
given by ~ ~ ) , so that the Hudson number must 
remain above one, thus excluding the field rotating solution of 
equations (41) and (42). The third root of equation (67) occurs 
wh<=*n 


t 



( 69 ) 
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Growing solutic^v^ require > 0, so that 


I 

- 

A 




When = 0 this collapses to Saffman's condition 1 < 

(note that 1 = sec 9). For the field reversing solu- 
tions of equations (41) and (42), n remains close to unity, so 
that the second term in equation (70) is less than 7.ero. Then a 
growing solution requires - (1 + Y ^ ( I + ys ^)) > 0, which 
requires mjj ^ ^ ^ If /S • Tluis the uniform field solution in 
the rotated field region cannot grow, and the field rotating 
solutions admitted by the momentum conservation integrals do not 


exist. 


The second class of solutions, described by equations (43) - 


(46) , involves an A1 fven number of order sec = 1 + ^ . Thus 
the first term in equation (70) is greater th<an zero, so that a 
growing solution requires 
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Thin con'iition must hoVl on either side of a nonuniform field 
— that iS/ for the values of » 3*^*^ V found on the side in 

question. To examine the situation in the rotated region, 


equation (4G) may be written in terms of m ^ as 





(72) 


where 
y3* and 
appear 


use has been 


o 

mu 1 t i pi ie<] 


made of the f 
nterchangeable 

by Y Using 


act that Y ^ ~ ^ 

to terms of order ^ 
equation (72) , the 


2 


and that 
when they 
condition 


(71) becomes 



(73) 


A second constraint is placed 
region is taken into account, 
well, while the definition ^ = 
expression : 


on ^ when the initial 
The inequality (71) holds 
mo“/(l +/? J) - 1 leads 


uniform 
here as 
to the 
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rif\ 




(74) 


Using equation (74) in inequality (71) leads to the condition 



(75) 


Comparison of i nequa I i t ier, (73) and (75) shows that growing solu- 
tions are oossihle in both the initial and rotated field regions 
provided 



(76) 


t 
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Inequality (76) shows that in warm plasma (^ >0) growing 
solutions are possible for both uniform field solutions admitted 
by the momentum conservation equations. However, the conditions 


of momentum conservation and asymptotic growth do not guaranty 
that two different uniform field solutions corresponding to the 
same initial conditions are actually connected. In the warm 
plasma case the existence of such a connecting solution can be 
determined only by numerical • integration, since the analytic 
integrals o^ the system are not sufficient to specify as a 
function of as is possible with cold plasma. 

Extensive numerical integration of the equations for systems 
of three to fifteen beams has confirmed the limitations on growth 


sugqeste«l by the asymptotic analysis. The resulting B , B 

^ y 

hodoqrams resemble those of Raffman even in the extended modified 


Alfven number range given by inequality (76). Although the 
parameters "D o' been varied extensively within 

aporopriate ranges, no solutions have been found which take one 

• 

uniform field solution into another with different field angle. 
Instead, true soliton solutions occur which are qualitatively no 
different from those of Saffman and which take the initial 


uniform fieirl solution back to itself. At the midpoint of the 
soliton, the y c«)mponent of the magnetic field vanishes but its 
derivative is always a maximum. Furthermore, the x component at 
this point does not have the appropriate rotated value given by 
equation (41). Thus the warm plasma theory appears incapable of 
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producing the rlonired ‘Boliitions. 

CONCLUSION 

The inclusion of thermal effects through use of a multibeam 
distribution function enriches and complicates the description of 
large amplitude disturbances. The extra degree of freedom 
afforded by the inclusion of the temperature-like variable Y 
a1 lows two classes of field rotating solutions to satisfy the 
momentum conservation integrals of the system. However, only 
those solutions associated with small changes in field angle 
appear capable of asymptotic growth on lx>th sides of the 
disturbance, and even these have not been found after extensive 
numerical investigations. It is unclear whether the apparent 
nonexistence of field reversing disturbances is a peculiarity of 
the particular and somewhat artificial model used or rather a 
reflection of the impossibility of such disturbances in colli- 
sion! ess pi asma . 
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Figure Caption 


Figure 1: A particle in a helical trajectory about the 

magnetic field. All space is assumed filled with such particles, 
whose orbits are in phase along the z axis, not along the field 
direction (z* 


axis) . 
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